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ABSTRACT 



We present results on Mrk 573 obtained as part of the CHandra survey of Ex- 
tended Emission-line Regions in nearby Seyfert galaxies (CHEERS). Previous studies 
showed that this source features a biconical emission in the soft X-ray band closely 
related with the Narrow Line Region as mapped by the [O m] emission line and the 
radio emission, though on a smaller scale; we investigate the properties of soft X-ray 
emission from this source with new deep Chandra observations. Making use of the 
^- subpixel resolution of the Chandra/ ACIS image and PSF-deconvolution, we resolve 

and study substructures in each ionizing cone. The two cone spectra are fitted with 
photoionization model, showing a mildly photoionized phase diffused over the bicone. 



Thermal collisional gas at about ~ 1.1 keV and ~ 0.8 keV appears to be located be- 



-^ tween the nucleus and the "knots" resolved in radio observations, and between the 

"arcs" resolved in the optical images, respectively; this can be interpreted in terms of 

,_i shock interaction with the host galactic plane. The nucleus shows a significant flux 

j> decrease across the observations indicating variability of the AGN, with the nuclear 

region featuring higher ionization parameter with respect to the bicone region. The 
long exposure allows us to find extended emission up to ~ 7 kpc from the nucleus 
along the bicone axis. Significant emission is also detected in the direction perpendic- 
ular to the ionizing cones, disagreeing with the fully obscuring torus prescribed in the 
AGN unified model, and suggesting instead the presence of a clumpy structure. 

Subject headings: galaxies: active - galaxies: individual (Mrk 573) - galaxies: jets - 
galaxies: Seyfert - X-rays: galaxies 
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1. INTRODUCTION 

The interaction between active galactic nuclei (AGN) and their host galaxy is signaled by 
the narrow line region and extended narrow line region (NLR, ENLR). This [O m] emitting re- 
gion, extending on kpc scale and observed in many Seyfert galaxies, is widely interpreted as gas 
photoionized by the AGN (e.g., Schmitt et al. 2003). 

A key tool to understand this interaction is provided by the study of soft X-ray emission 
from these sources. The overall structure of the X-ray emission appears to be morphologically 
correlated with the ENLR, suggesting a common physical origin (Young et al. 2001; Iwasawa et 
al. 2003; Bianchi et al. 2006; Wang et al. 201 lc). 

In the unified AGN model (Antonucci 1993) the soft X-ray spectra from Seyfert 2 galaxies 
are affected by emission from the circumnuclear medium, illuminated by the nuclear continuum. 
This is supported by high resolution spatial/spectral observations with both Chandra and XMM- 
Newton showing that soft X-ray emission arise from the circumnuclear medium photoionized by 
the nuclear continuum. Moreover, X-ray grating spectroscopy agrees with this picture, showing 
that the spectra can be interpreted as a blending of emission lines, with little contribution from 
collisionally ionized plasma (see e.g. Ogle et al. 2000; Brinkman et al. 2002; Kinkhabwala et al. 
2002; Guainazzi & Bianchi 2007). 

Mrk 573 (z = 0.017, Ruiz et al. 2005) is one of the brightest Seyfert galaxies observed by HST, 
with a ENLR extending to ~ 9" corresponding to a projected size of about 3 kpc 1 , where photoion- 
zation by the central AGN appears to be the dominant process (Ferruit et al. 1999; Schlesinger et 
al. 2009). Mrk 573 is also also associated with a triple radio source, with a central core and two 
side knots (Ulvestad & Wilson 1984). 

The ENLR of Mrk 573 has been modeled by Fischer et al. (2010) making use of HST STIS 
long- slit spectra; the kinematic model developed by these authors features an ionizing bicone with 
the NW cone tilted toward the observer line of sight, and with a half opening angle of 53° which 
is close to the 60° expected in the unifying structure for the inner quasar regions presented by 
Elvis (2000). The observed half opening angle of ~ 30° results from the intersection of the bicone 
with the galaxy plane. The authors conclude that the circumnuclear emission is mainly due to the 
intersection of the ionizing bicone with the galactic disk, as supported by the observed emission 
arcs that coincide with outer dust lanes. 

The X-ray emission from Mrk 573 has been previously studied by several authors. A study 



'in the following, we adopt the standard flat cosmology with Qa = 0.73 and Hq = 70 km s ' Mpc ' (Komatsu et 
al. 2011). 
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of EPIC/XMM-Newton data by Guainazzi et al. (2005) showed for this source a strong Fe Ka line 
with an equivalent width larger than 1 keV, yielding to a classification as a Compton-thick object; 
high-resolution spectral analysis of RGS/XMM-Newton data by Guainazzi & Bianchi (2007) led 
to an interpretation of the X-ray emission from this source in term of gas photoionized by the cen- 
tral AGN. Recently, two papers concerning the study of soft X-ray emission from Mrk 573 with 
ACIS/ 'Chandra and RGS/XMM-Newton data have been published (Bianchi et al. 2010; Gonzalez- 
Martin et al. 2010); they both agree in interpreting this emission as dominated by two photoionized 
phases, with the contribution of a collisional phase. Making use of the archival Chandra observa- 
tion available in 2010, these authors found a clear correlations between the diffuse X-ray emission 
from Mrk 573, its ENLR as mapped by HST data, and the triple radio source shown in the VLA 
image, suggesting a physical interplay among these three components. 

In this paper we present a detailed imaging and spectroscopic study of new Mrk 573 AClS/Chandra 
data obtained as part of the CHandra survey of Extended Emission-line Regions in nearby Seyfert 
galaxies (CHEERS, Wang 2010). These new data, adding to the archival ones, reach an observa- 
tion time of ~ 110 ks. With this deep observation the unmatched Chandra spatial resolution allows 
us to explore different regions of the Mrk 573 ENLR looking for a more comprehensive description 
of the underlying physical processes. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. AClS/Chandra Data 

Mrk 573 was observed by Chandra on 2006 November 18 for an exposure time of 40 ks (Obs. 
ID 07745, PI: Bianchi), and on 2010 September 16, 17 and 19, for an exposure time of 10, 53 and 
17 ks respectively (Obs. IDs 12294, 13124, 13125); the last three observations were performed as 
part of the CHEERS survey. Level 2 event data were retrieved from the Chandra Data Archive 2 
and reduced with the Chandra Interactive Analysis of Observations (CIAO; Fruscione et al. 2006) 
4.3 software and the Chandra Calibration Data Base (caldb) 4.4.3, adopting standard procedures. 

Time intervals of background flares exceeding 3<x of the quiet level were excluded using the 
lc_sigma_clip task in source free regions of each observation; net exposure times are reported in 
Table 2. The nucleus has no significant pile up, as measured by the CIAO pileup_map tool. 

We produced a merged image of the four observations, to take advantage of the longer expo- 
sure time and identify fainter signatures; to this end we used the wavdetect task to identify point 

2 http://cda. harvard.edu/chaser 
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sources in each observation, then we used the reproject_aspect task to modify the aspect solution 
minimizing position differences between the sources found, and finally merged the images with 

the MERGE_ALL SCHpt. 



2.2. Optical and Radio Data 

To compare the soft X-ray emission with optical and radio structures, we retrieved HST 
and VLA images of Mrk 573. In particular, we retrieved from the Hubble Legacy Archive 3 the 
//5T/WFPC2 narrow band image with the FR533N filter tracking the [O m]A5007 emission line, 
obtained as part of the GO program 6332 (PI: Wilson) on 1995 November 12 (Falcke et al. 1998; 
Schmitt et al. 2003); this image was aligned to the X-ray images by matching the brightest pixel of 
the nuclear source in both images. The [O m] image (Figure 1, upper-right panel) shows the bright 
nucleus with two bright spots at ~ 1" and ~ 0.7" in the NW and SE directions, respectively; we 
also see two pairs of arcs in the two cones, the brighter inner arcs at ~ 2" and ~ 1 .5" and the fainter 
outer arcs at ~ 4" and ~ 3" from the nucleus in the NW and SE directions, respectively (see also 
Ferruit et al. 1999; Quillen et al. 1999). 

We retrieved from the NASA/IPAC Extragalactic Database (NED 4 ) the VLA image of the 
6 cm (5 GHz) observation performed in 1985 March 3 in A array configuration; this image clearly 
shows a central nucleus with a size ~ 1" and two side lobes of similar size aligned with the NLR 
(Figure 1, contours in lower-right panel, see also Ulvestad & Wilson 1984; Falcke et al. 1998). 



3. CHANDRA IMAGE ANALYSIS 

Imaging analysis was performed with the subpixel event repositioning (SER) procedure (Li 
et al. 2003) and without pixel randomization, to take advantage of the telescope dithering to allow 
subpixel binning of the images, using a pixel size smaller than the native one 0.492" of the Chan- 
dra/AClS detector (see, e.g., Harris et al. 2004; Siemiginowska et al. 2007; Perlman et al. 2010; 
Wang etal. 2011a). 



3 http://hla. stsci.edu/hlaview.html 
4 http://ned. ipac.caltech.edu/ 



3.1. Comparison with Chandra PSF 

We performed Point Spread Function (PSF) simulations with Chandra Ray Tracer (ChaRT 5 , 
Carter et al. 2003) taking into account the source spectrum, exposure time and off-axis angle; while 
the hard band images (2-10 keV) show a point-like structure consistent with the simulated PSF, 
the soft band images (0.3 - 2 keV) show extended morphology that closely resembles that of the 
NLR, as mapped by the [O m] emission (see Figure 1). 

In Figure 2 we show radial profiles of the soft X-ray emission for the NW and SE direction 
compared with the simulated PSF in the same band; while the nuclear emission in the inner ~ 1" 
is comparable with the PSF, the diffuse emission extends in both directions up to ~ 12". 



3.2. Image Deconvolution 

We applied to the merged images of the soft emission two different PSF-deconvolution al- 
gorithms, namely the Richardson-Lucy (R-L) (Richardson 1972; Lucy 1974) and the Expectation 
through Markov Chain Monte Carlo (EMC2) (Esch et al. 2004; Karovska et al. 2005, 2007); while 
both methods show a similar extended morphology in the deconvolved images, the R-L algorithm 
yields a more grainy image reconstruction with respect to EMC2; the latter results are presented in 
Figure 1 . In particular, in the upper panels we show a comparison between the reconstructed ACIS 
image with the subpixel binning (1/8 of the native pixel size, upper-left panel) and the [O m] image 
(upper-right panel): the soft X-ray extended emission clearly shows structures in correspondence 
with the optical optical arcs (Hollis et al. 1997); this correspondence is more clearly shown in the 
lower panels of the same figure, where [O m] and radio contours are overlaid on the soft X-ray 
emission (lower-left and lower-right panel, respectively). While the radio emission appears more 
compact than the soft X-ray emission, the latter shows a striking coincidence with the optical fea- 
tures and an interesting interplay between radio, optical and X-ray emissions, as already suggested 
by Bianchi et al. (2010) and Gonzalez-Martin et al. (2010) in their analysis of Chandra obs. 07745. 
In particular we see soft X-ray structures lying in coincidence or just in front of the inner [O in] 
arcs, with the latter wrapping around the outer radio lobes. 



5 http ://cxc .harvard.edu/chart/ 



3.3. Diffuse emission 

The deep exposure reached with CHEERS observations (~ 1 10 ks) allows us to explore fainter 
features of the diffuse X-ray emission. To enhance these, we applied the adaptive smoothing pro- 
cedure to the merged soft X-ray image using the csmooth tool (Ebeling et al. 2006), with minimum 
and maximum significance S/N levels of 2.5 and 3.5, respectively; this smoothing procedure allows 
us to enhance fainter, extended features of the diffuse emission. The smoothed image is presented 
in Figure 3; the bicone emission appears to extend to ~ 6 kpc in the NW direction and ~ 7 kpc in 
the SE direction, as also confirmed by the extended radial profiles presented in Figure 4. 

We then focused on the region perpendicular to the bicone that is supposed to be shielded by 
a dusty torus from seeing the central continuum in the AGN unified model; the region we studied 
is defined by a bicone with P.A. = 34°, a half opening angle of 55°, and inner and outer radii of 
1.5" and 12", as shown in Figure 5. We found a significant (SNR « 22) soft X-ray emission in this 
region, extending to ~ 9" as shown in Figure 6; we will discuss the emission from this region in 
detail in Sect. 5. 

The nuclear region has a luminosity in the 0.3 -2 keV band of 24.16^oio x ^ 4 ° er S s_1 > wn il e 
the bi-cone region has a luminosity of 3.26^;[2 X 10 40 erg s" 1 and 2.97^;}f x 10 40 erg s" 1 in the NW 
and in the SE direction, respectively. The cross-cone region, on the other hand, has a luminosity of 
l^gxlO^ergs- 1 . 



4. SPECTRAL ANALYSIS 

In order to study the X-ray properties of the diffuse emission we extracted spectra with the 
CIAO specextract task from three different regions centered at RA 01 :43:57.78, DEC +02:20:59.32: 
the NW and SE cones region, both defined as cones with an inner radius of 1.5" and a half opening 
angle of 30° , centered at P.A. = -56° and 1 24° respectively; the nuclear region, defined as a circular 
region with a 1" radius (see Figure 5). Based on the radial profiles shown in Figure 4 we decided 
to extract cones spectra out to 18", where the emission of the cones reaches the background level. 
For spectra extracted from the nuclear region we applied the point-source aperture correction to 

SPECEXTRACT task. 

To make use of the;^ 2 fit statistic we binned the spectra to obtain a minimum of 20 counts per 
bin using the specextract task; in the following, errors correspond to the 1-cr confidence level for 
one interesting parameter (A^ 2 = 1). In all the spectral fits we included photo-electric absorption 
by the Galactic column density along the line of sight N H = 2.52 x 10 20 cm" 2 (Kalberla et al. 2005). 
We also tried to evaluate intrinsic absorption, however the data do not show a significant intrinsic 
column density. 



-7 



4.1. Cone spectra 

The soft X-ray emission from Seyfert galaxies can be effectively described in terms of several 
emission lines with a small power-law continuum contribution. We therefore performed the spec- 
tral analysis of emission of Mrk 573 with the XSPEC software (ver. 12.7.0 6 , Arnaud 1996) fitting 
the soft (0.3 - 2 keV) spectra of different observations in the cones region using a phenomenolog- 
ical model constituted by a power-law photon index fixed to Y = 1.8 (a typical value for Seyfert 
galaxies, see Bianchi et al. 2009), plus several red-shifted emission lines adopting as a reference 
the lines measured in the 2006 November 18 observation (Obs. ID 07745) by Gonzalez-Martin et 
al. (2010) in RGS/XMM-Newton high resolution spectra; due to low counts in these regions, we 
were able to obtain well constrained fits only for longest observations, namely Obs. IDs 07745 
and 13124. To obtain better statistics, we also performed simultaneous fits on data from CHEERS 
observations (Obs. IDs 12294, 13124 and 13125) as well as for all observations (CHEERS obs. + 
Obs. ID 07745). To avoid contamination from the nuclear continuum we evaluated the contribu- 
tion from the PSF wings in the cone regions as ~ 0.8%, and subtracted this contribution from the 
cones spectra. Fit results are in general agreement with previous results by Bianchi et al. (2010) 
and Gonzalez-Martin et al. (2010), and are presented in the Appendix in Table 4 and Figure 12. 

The strongest lines of the cones spectra are C vi, O vn triplet, O vm, O vn RRC and Ne ix 
triplet. Most lines have similar flux in both cones, except for O vm RRC which appears stronger 
in the NW cone, as already reported by Gonzalez-Martin et al. (2010), and the Ne ix triplet which 
instead is stronger in the SE cone. 

The high resolution RGS/XMM-Newton spectroscopic analyses of Bianchi et al. (2010) and 
Gonzalez-Martin et al. (2010) suggest that the soft X-ray spectrum of Mrk 573 originates from 
photoionized and photoexcited plasma with a contribution from collisionally excited plasma. In 
order to evaluate the contribution of different plasma phases to the total emission we performed 
spectral fitting of the two cones spectra separately, making use of self-consistent photoionization 
models. For this purpose we produced xspec grid models with the Cloudy 7 c08.01 package, de- 
scribed by Ferland et al. (1998). We assumed the ionization source to be an AGN continuum (with 
a big bump temperature T = 10 6 K, a X-ray to UV ratio a ox = -1.30 and an X-ray power-law 
component of spectral energy index a = -0.8) illuminating a cloud with plane-parallel geometry 
and constant electron density n e = 10 5 cm -3 . Note that the fits are expected to be quite insen- 
sitive to n e in this density regime (Porquet & Dubau 2000). The grid of models so obtained are 
parametrized in terms of the ionization parameter U (varying in the range log U = [-3.00 : 2.00] 



6 http://heasarc. nasa.gov/xanadu/xspec/ 
7 http://www.nublado.org/ 



in steps of 0.25) and the hydrogen column density N H (expressed in cm" 2 varying in the range 
logN H = [19.0 : 23.5] in steps of 0.1), taking into account only the reflected spectrum from the 
illuminated face of the cloud (Bianchi et al. 2010; Marinucci et al. 201 1). 

The results of the fitting procedure are shown in the Appendix in Table 5 and Figure 13. Poor 
fits are obtained with a single phase ionized plasma; the best fit models require two photoionized 
phases. Due to low counts, the hydrogen column densities of both phases are unconstrained, so we 
fixed them to their best fit value of log Nh = 20. 

The simultaneous fit of all observations give for the NW cone log Ui = 0.9 + 0.3 and log U2 = 
-0.5 + 0.3 for the higher and lower ionization phase, respectively; for the SE cone, on the other 
hand, we have log U\ = 0.7 + 0.2 and log U 2 = -0.8 + 0.2. Comparable values are obtained for the 
simultaneous fit of all the CHEERS observation, as well as for obs. ID 13124 and 07745, even if 
for the latter the ionization parameters log U2 have to be fixed to their best fit values -0.5 and -1 
for the NW and SE cone, respectively. In all observations the two gas phases show similar fluxes 
F = (0.3 + 0.1) x 10" 13 erg cm" 2 s" 1 in both cones. 

In almost all observations the best fit models require a thermal component (included with an 
apec xspec model). The exception is observation 07745 on the SE cone, but this is probably due to 
the low statistics, this being the region with the lowest counts. All other fits require this thermal 
component with temperature kT = 0.8 ±0.1 keV in the NW cone and kT = 1.0 + 0.2 keV in the 
SE cone. 

We can therefore conclude that a two-phases photoionized plasma is diffused over the two 
cones, with the presence of a thermal component. 



4.2. SPATIAL DISTRIBUTION OF SPECTRAL FEATURES 

Summarizing the results obtained in the spectral fits presented in Sect. 4.1 the two cones are 
characterized by mildly photoionized plasma (log U\ = 0.9 + 0.3 in the NW cone and log U\ = 
0.7 + 0.2 in the SE cone) and the presence of collisionally ionized plasma (at a temperature kT = 
0.8 + 0.1 keV in the NW cone and Jfcr = 1.0 ± 0.2 keV in the SE cone). 

As recalled in Sect. 3.2 CHEERS observations of Mrk 573 reveal a detailed structure of 
the diffuse emission that is morphologically coincident with the ENLR, as mapped by [O m] and 
suggest an interplay between radio ejecta, the optical and the soft X-ray emissions (see Figure 1); 
the inner optical arcs appear wrapped around the radio knots, as resulting from shock interactions 
between the radio jets and the galactic plane, and the X-ray emission is enhanced in front of these 
knots, suggesting the presence of collisionally ionized gas, a picture analogous to what is observed 



in NGC 4151 (Wang et al. 2011b). 

To further investigate these features, we produced maps of the X-ray emission in three energy 
bands corresponding to some of strongest emission lines, namely the blended O vn triplet (0.53 
- 0.63 keV), O vm Lya (0.63-0.68 keV) and Ne ix triplet (0.90 - 0.95 keV), and present them in 
Figure 7. 

As already observed by Bianchi et al. (2010) the O vn emission is slightly more extended in 
the SE cone with respect to the O vm Lya, while in the NW cone we observe an opposite behavior. 

The Ne ix map shows spots of enhanced emission in front of the radio lobes; these spots are 
also present in the ratio map of line emission shown in Figure 8, where we present ratios of Ne 
ix to O vn and to O vm, tracing the higher ionization gas; these maps show a nearly symmetrical 
structure along the two cones, characterized by regions of high ionization between the nucleus and 
the inner optical arcs (regions 1 and 2) and regions of increased ionization lying in front of the 
radio knots at the interface between the inner and outer optical arcs (regions 3 and 4). 

These regions of increased ionization are also shown in Figure 9, where we plot variations of 
the ratio of [O m] to 0.5 - 2 keV flux along with the distance from the nucleus, where an average 
conversion factor of 8.7 x 10~ 12 erg cm" 2 s" 1 counts" 1 s, obtained by the various spectral fits we 
performed, is used to convert net counts in the 0.5 - 2 keV band into flux. The ratio of [O in] to 
soft X-ray flux for a single photoionized medium is expected to have an approximatively power 
law dependence on the radius, depending on the radial density profile (Bianchi et al. 2006). The 
trends shown in Figure 9 can be roughly sketched with the inner ~ 1" region corresponding to 
the nuclear highly photoionized phase with log U = 1, and two outer regions corresponding to 
the mildly photoionized cones with log U = 0.8 in the NW direction and log U = 0.6 in the SE 
direction, respectively. 

As the radio ejecta interact with the galactic plane, the regions of increased ionization shown 
in Figure 8 can be interpreted as two different shock regions in which the radio ejecta interact with 
different regions of the ISM, with decreasing velocities as the jet slows down moving away from the 
nucleus; in particular, regions 3 and 4 feature significant dips in the [O m] to soft X-ray flux ratio 
at the interface between the optical arcs (see also Gonzalez-Martin et al. 2010), characterized by 
enhanced X-ray emission due to shock heating in addition to photoionization, as already observed 
in NGC 4151(Wang et al. 2011b). 

In order to investigate the nature of the different plasma phases in these regions of enhanced 
ionization (shown in Figure 8) we extracted here soft X-ray spectra and fitted them adopting the 
photionization model described above, taking into account the nuclear contribution to these re- 
gions. In order to have enough counts and allow acceptable fits we extracted spectra from all 
observations and combined them with the specextract tool. Due to the low counts in regions 3 and 



-10 



4, we combined their spectra. Fits results are presented in Table 3 and Figure 10. 

In regions 1 and 2 the best fit model requires a highly ionized phase with log U\ = 0.9 ±0.1 
and Nm = 21.5 ± 0.6, and a low ionized phase with log U<i = -1.2 ± 0.6 in region 1 and log U2 = 
-0.8 ± 0.3 in region 2; due to low statistics the hydrogen column densities for the low ionization 
phase were fixed to their best fit values of logN H2 = 20 and \ogN m = 21 in region 1 and 2, 
respectively. The spectrum extracted from regions 3+4 is best fitted by a highly ionized phase with 
log U\ = 1.3 ± 0.3 and a low ionized phase with log U% = -0.9 ± 0.2; in this case both hydrogen 
column densities were fixed to their best fit value log Nh = 20. For all these regions of enhanced 
ionization a thermal component is statistically required; in particular, region 1 needs a temperature 
kT = 1.1 ± 0.2 keV, while in region 2 we have kT = 1.3 ± 0.2 keV. The outer regions 3+4 feature 
a lower temperature plasma with kT = 0.8 ±0.1 keV, which we interpret as the origin of the 
enhanced X-ray emission. 

To further check the location of the collisionally ionized phases we performed an analogous 
analysis on the [O m] outer arcs. A good fit is obtained with two photoionized phases, while the 
inclusion of a thermal component is not statistically required, as the fitting procedure only yields 
an upper limit on its flux < 0.02 x 10 erg cm" 2 s" 1 (see Table 3 and Figure 10). This is not 
compatible with the fluxes found in the regions of enhanced ionization (see Table 3). We conclude 
that the thermal component found by Bianchi et al. (2010) and Gonzalez-Martin et al. (2010) is 
actually composed of two physically separated plasmas, with a higher temperature phase lying 
between the nucleus and the inner optical arcs, and a lower temperature phase located in front of 
the radio knots, at the interface between the inner and outer optical arcs. 



5. DISCUSSION 

CHEERS observations of Mrk 573 show soft X-ray emission characterized by a bright nu- 
cleus and a diffuse component extending along the ionizing bicone up to ~ 7 kpc, doubling the 
size reported in previous works (Gonzalez-Martin et al. 2010). The spectral analysis of the diffuse 
emission points to the presence of a thermal component that can be interpreted as collisionally ion- 
ized gas related to radio ejecta; moreover, significant extended emission up to ~ 3 kpc is observed 
in the cross-cone direction where the obscuring torus is expected to shield the nuclear continuum 
according to the AGN unified model. In this section we will discuss these two aspects of the diffuse 
emission, as well as the spectral properties of the nuclear emission. 



-11 



5.1. Energetic content of the collisionally ionized gas 

Previous work on Mrk 573 (Gonzalez-Martin et al. 2010; Bianchi et al. 2010) first suggested 
that the soft X-ray diffuse emission of Mrk 573 can be interpreted as originating from mildly 
photoionized plasma with the presence of a collisionally ionized phase. With new, deep CHEERS 
observations we are able to conclude than the collisionally ionized phase (Sections 5.2 and 4.1) 
is composed of two physically separated plasmas, with a hotter plasma between the nucleus and 
the inner optical arcs, and a lower temperature plasma in front of the radio knots, at the interface 
between the inner and outer optical arcs. 

The collisionally ionized phases can be interpreted as the result of shock interaction between 
the radio knots and the galactic plane, because this is the expected location of the interaction 
between the ionizing cones and the galaxy plane (see Figure 1). 

The electron density and thermal pressure of the collisionally ionized gas can be estimated 
through the normalization of the apec component, that is, its emission measure 



EM = 10 



-14 n e n H Vri 
4tt[D a (1 + z)f 



where n e and n H are the average electron and hydrogen density, V is the emitting region volume, 
rj is the filling factor and D A is the angular diameter distance to the source. The gas parameters 
were obtained assuming n e = 1 .2 n H , a filling factor r\ « 1 and a spherical shape for the emitting 
regions. All parameters are listed in Table 1. As shown, the electron density, the thermal pressure 
and the temperature decrease from the inner (1 or 2) to the outer shock regions (3+4). Considering 
an average electron density over all regions of increased ionization (n e ) =0.18 cm -3 and a gas 
mass-averaged temperature (T) = 0.88 keV, we can evaluate an average thermal pressure (p th ) = 
5.7 x 10~ 10 dyne cm" 2 . If the collisionally ionized gas is assumed to be diffuse over the entire cone 
region a thermal pressure 1.5 x 10~ 9 dyne cm -2 is found, in agreement with the result found by 
Bianchi et al. (2010). 

To compare the thermal pressure with the pressure of radio jets we assume energy equiparti- 
tion between particles and magnetic field in the jet, which provide a lower limit to the jet pressure 
(see, e. g., Hardcastle et al. 2004). We model the two radio knots as spheres with a radius ~ 200 pc 
filled with electrons having a power law energy distribution with an index r = 2.4 (the average of 
the indices measured by Falcke et al. 1998), and Lorentz factors 2 < y < 10 5 . Assuming that the 6 
cm luminosities reported by Bianchi et al. (2010) are due to synchrotron emission from this elec- 
tron population we obtain values of the magnetic field of 6.6 x 10~ 9 G and 8.0 x 10~ 9 G for the NW 
and SE lobe, respectively, and jet pressures due to electrons energy density of 3 x 10 -11 dyne cm -2 
and 4 x 10~ n dyne cm -2 for the NW and SE lobe, respectively. These are about 10 times smaller 
than the thermal pressure in regions 3+4, ~ 60 times smaller than the thermal pressure in regions 
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1 and 2, and ~ 20 times smaller than the average thermal pressure evaluated in all regions of in- 
creased ionization. Therefore the jet pressure due to electrons is not sufficient to sustain shock 
heating of the hot gas (Bianchi et al. 2010). 

However, if we consider jet pressures due to the total energy density, 7/3 t/g (where Ub is the 
magnetic field energy density), we obtain 1.5 X 10~ 10 dyne cm" 2 and 2.2 x 10~ 10 dyne cm" 2 for the 
NW and SE lobe, comparable with both the average thermal pressures (p t h) and the thermal pres- 
sure of the lower temperature gas found at the interface between the optical arcs (region 3+4). This 
suggests that shocks of radio ejecta with ISM are likely to be responsible for the gas heating we 
have located in front of the radio lobes. Nevertheless, due to large uncertainties in these estimates, 
a strong conclusion cannot be drawn. 

The thermal energy associated with the collisionally ionized gas in the regions discussed 
above is E thw = 7.9 x 10 53 erg, E th{2 ) = 1.0 x 10 54 erg E t h(3+4) = 2.7 x 10 54 erg for regions 
1, 2 and 3+4, respectively; the corresponding cooling times t coo i = E t h/L c (where L c is the lu- 
minosity of the thermal component) are ~ 10 7 yr. In these regions the local sound speed c s 
is comparable with the shock velocity estimated assuming strong shock conditions with v s h ~ 
100 (£770.013 keV) 1/2 km s" 1 (Raga et al. 2002), resulting in crossing times of the radio knots 
tcross ~ 10 5 yr, about 30 times shorter than the cooling time. Thus no extra heating sources are 
required. 

A lower limit on the fraction of the kinematic power converted to heat in the ISM can be 
estimated as L K > E th /t cross , to yield L K \ > 6.9 x 10 40 erg s _1 , L K2 > 1.0 x 10 41 erg s _1 and 
Lk(3+4) ^ 2.0 x 10 41 erg s" 1 for regions 1, 2 and 3+4, respectively. This can be compared with 
the jet total power, that can be evaluated using the Pj et - P radio relation in Cavagnolo et al. (2010) 
and the 1.4 GHz radio luminosity vL v = 4.1 x 10 37 erg s" 1 (Condon et al. 1998), to obtain P jet « 
1.3 x 10 42 erg s -1 . This yield that > 5% of the jet power is deposited into region 1, > 7% is 
deposited in the ISM into region 2 and > 15% is deposited into region 3+4, so that a total fraction 
of > 28% of the jet power is deposited in the ISM, a fraction much larger than the 0. 1% estimated 
in NGC 4151 (Wang et al. 201 lc). Considering the ANG bolometric luminosity evaluated from the 
2-10 keV emission (see Sect. 5.2) L bol = (4.4 + 7.7) x 10 44 erg s" 1 we can conclude that a fraction 
L K /L bo [ > 0.05% of the available accretion power is thermally coupled to the host ISM. This is 
substantially lower than the fraction of the accretion power of the black hole (~ 5% - 100% of Lboi 
) that the majority of quasar feedback models require to be deposited into the host ISM (Mathur et 
al. 2009), but is somewhat compatible with a two-stage feedback model proposed by Hopkins & 
Elvis (2010), which requires only 0.5% of the radiated energy to drive the mass outflow. 
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5.2. Nuclear properties 

In this section we focus on the nuclear emission of Mrk 573. The region used for the spectral 
extraction, shown in Figure 5, is a circle centered in RA 01:43:57.78, DEC +02:20:59.32 with a 
1" radius. 

To evaluate the emission lines fluxes we adopted the same phenomenological model used in 
Sect. 4.1; these results also confirm results by Bianchi et al. (2010) and Gonzalez-Martin et al. 
(2010), and are presented in the Appendix in Table 6 and Figure 14. 

The strongest lines of the nuclear spectra are C v Hey, C vi Ly/3, N vn Lya, O vn triplet, 
O vm Lyar, O vn RRC, Fe xvn and Ne ix triplet. Significant variations are observed between 
observations. In particular, observed continuum flux in the 0.3 - 2 keV band f sx is higher by a 
factor 1.3 x 10" 13 erg cnr 2 s" 1 (4cr) in OBS. 07745 [fsx = (5.4 + 0.3) x 10" 13 erg cnr 2 s" 1 ] than in 
the 4 years later CHEERS observations [f sx = (4.1 + 0.1) x 10~ 13 erg cm" 2 s" 1 ]. These variations 
are reflected in the poor fit obtained with merged observation (x 1 = 1.63). To show the observed 
nuclear flux variations we plot in Fig. 1 1 fluxes of the different spectral components used in our 
phenomenological model. We clearly see that the nuclear flux decrease is marked by a significant 
decrease in several emission lines: O vn triplet, O vm Lya, O vn RCC and Ne x Lya, at 3cr. 

We then fit the nuclear spectra with the photoionization model described in Sect. 4.1. Since 
the nuclear emission significantly extends beyond the soft band, we performed nuclear spectral fits 
over the full 0.3 - 10 keV energy range, adopting, in addition to the photoionization model: a pure 
neutral reflection component with the pexrav xspec model with the photon index fixed to 1.8, cos i 
fixed to 0.45 and Si Kor, S Ka and Fe Ka lines (see Bianchi et al. 2010). The results of the fitting 
procedures are presented in the Appendix in Table 7 and Figure 15. 

As for the cones, poor fits are obtained with a single phase ionized plasma; the best fit models 
require two photoionized phases, with higher and lower values of the ionization parameter. 

We also tired adding a thermal component, included with the apec xspec model, looking for 
the thermally ionized plasma found by Gonzalez-Martin et al. 2010 and Bianchi et al. 2010. The 
thermal component, however, is not statistically required by the best fit models and does not pro- 
vide better x 2 values, yielding only upper limits for the thermal component flux F c . We find 
for CHEERS obs. Fc(o.3-io) < 0.06 x 10 _13 erg cnr 2 s" 1 , while for all merged obs. we find 
Fc (0.3-10) < 0.08 x 10~ 13 erg cm" 2 s" 1 . These fluxes are comparable with those of the hot colli- 
sionally ionized plasma highlighted in the ionization maps presented in Figure 8, and due to the 
strong nuclear emission it is not possible to constrain its parameters when extracting spectra in 
the central 1" region, so we included it in the fitting model as a constant component taking into 
account the contribution of regions 1 and 2 to the nuclear region (see Figure 8). 
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The higher ionization phase features values of log U\ = 1.1+0.1 and \ogN m = 20.6 + 0.2 
for the ionizing parameter and the hydrogen column density, respectively, both in obs. 07745 and 
in CHEERS observations; the best fit values for the less ionized phase are log U 2 = -0.75 + 0.2 in 
all observations, and log N m = 20.2 ± 0.4 and log N m = 21.0 ± 0.6 in obs. 07745 and in CHEERS 
observations, respectively. 

The photoionization model provides soft X-ray fluxes consistent with those obtained with the 
phenomenological model described above: for obs. 07745 fsx = (5.8 + 0.2) x 10~ 13 erg cm' 2 s -1 , 
and for the CHEERS obs f sx = (4.1 + 0.1) x 10~ 13 erg cm' 2 s' 1 . The flux in the 2 -10 keV 
band f HX is f HX = (3.8 + 0.3) x 10~ 13 erg cm" 2 s" 1 for obs. 07745 and f HX = (3.5 + 0.3) x 
10~ 13 erg cm" 2 s" 1 for CHEERS obs., which are consistent within uncertainties. The ionization 
parameter and the hydrogen column density for the two sets of observations are also consistent 
within the uncertainties. These results are confirmed by spectral fits to a restricted nuclear region 
(namely, the inner 0.5" radius of the nuclear region, shown in Figure 8) that does not include the 
regions 1-2 of increased ionization, and whose results are presented in Table 8 in the Appendix. 
These fits also show a decrease of f sx and a constant f HX across the observations, together with 
comparable parameters of the photoionization model. 

The simplest interpretation for the observed variations of the soft X-ray flux is for them to 
be related to intrinsic source variability, with the soft X-ray emission originating close to the line 
of sight. In this case the hard emission, constituted by reflected nuclear emission as expected 
for a Compton-thick source, must originate farther away from the line of sight. A possible site 
for the reflected emission to arise is the inner walls of the dusty torus (Ghisellini et al. 1994; 
Iwasawa et al. 1997; Matt et al. 1997; Bianchi et al. 2001) or warm absorber (see Sect. 5.3) lying 
at r > 4 ly = 1.2 pc from the nucleus, based on the time lag between observations. From f HX we 
evaluate an observed luminosity L 2 -w = (2.4 + 0.1) x 10 41 erg s" 1 , while the intrinsic luminosity 
can be recovered correcting this value by a factor 60 (Panessa et al. 2006) or 100 (Cappi et al. 
2006) due to the Compton-thick nature of Mrk 573 (with a column density Nh > 1.6 x 10 cm -2 , 
Guainazzi et al. 2005), to find L 2 -w = (1.5 4- 2.5) x 10 43 erg s -1 . Then we can evaluate the AGN 
bolometric luminosity using a bolometric correction of 30 (Risaliti & Elvis 2004), to find L ho i = 
(4.4 4- 7.7) x 10 44 erg s" 1 , consistent with L bo i = 5.2 x 10 44 erg s" 1 evaluated by Gonzalez-Martin 
et al. (2010) modeling the optical and infrared spectrum of Mrk 573 with a clumpy torus model, 
and slightly smaller than the L bo i = 7.9 x 10 44 4 1.3 x 10 46 erg s' 1 derived by Kraemer et al. (2009) 
using the correlation between the [Oiv] and L 2 - W . The bolometric luminosity can be used to set a 
lower limit on the torus inner radius r = 0.4L 45 1/2 r 150 o~ 26 pc, where L 45 is the bolometric AGN 
luminosity in units of 10 45 erg s _1 and r 150 o is the dust sublimation temperature in units of 1500 K 
(Nenkova et al. 2008). Using our estimate of the bolometric luminosity and setting r = r^ we 
obtain a sublimation temperature ~ 1000 K characteristic of silicate dust grains. 
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5.3. Emission from the cross cone region 

As discussed in Sect. 3.2 we found a significant emission in the cross cone direction. This 
emission is not expected in the presence of a continuous obscuring dusty torus, and may therefore 
be due to nuclear radiation leaking through a clumpy structure (Nenkova et al. 2008). Cross-cone 
emission was also found in NGC 4151 by Wang et al. (2011c), where the authors interpreted it 
as the result of leaking nuclear continuum through a filtering structure with an uncovered fraction 
~ 1%. 

Fischer et al. (2010) developed a kinematic model of Mrk 573 NLR, and concluded that the 
apparent bicone half opening angle of ~ 30° is likely due to the intersection with the galactic plane 
of a much wider ionizing bicone, with an intrinsic half opening angle ~ 53°, close to the 60° for 
the warm absorber structure presented by Elvis (2000). 

In order to evaluate the uncovered fraction of this filtering structure we fitted the spectrum 
of the soft emission in the cross-cone region shown in Figure 5 with a photoionization model, 
taking into account contamination from nuclear emission evaluated from the contribution of the 
PSF wings in the cross-cone regions to be ~ 2.3%. 

We used the same grid models produced with the Cloudy package used above, adopting the 
same model adopted for fitting the cone emission with two photoionized phases; the result of this 
fit is presented in Table 3 and Figure 10. 

The hydrogen column densities of both phases have been fixed to their best fit value of 
logN H = 20; the higher ionization component has a ionization factor log U\ = 0.8 ±0.1 simi- 
lar to the value obtained for the cone region, while for the lower ionization component we put an 
upper limit on the ionization factor log U<i < -2.4. The latter value is significantly lower than the 
one we obtain in the cone region, indicating that a smaller fraction of the nuclear continuum leaks 
in the cross-cone region. In order to estimate the uncovered fraction of the warm absorber we can 
use the ionization parameters obtained for the cone regions and for the cross-cone region. Since 
U ~ Q/r 2 , where Q is the number of nuclear ionizing photons per second illuminating the clouds, 
and assuming the same distance r for the cones and the cross-cone region from the nucleus, the 
ratio of the ionization parameters gives us the ratio of the nuclear ionizing photons seen by the dif- 
ferent regions. So, if we assume the cone regions to see the unshielded nuclear radiation, this ratio 
gives for the shielding structure an uncovered fraction < 3%, consistent with the fraction evaluated 
in NGC 4151. 
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6. CONCLUSIONS 

We have presented results from new AClS/Chandra observations of Seyfert 2 galaxy Mrk 573 
obtained as part of CHEERS survey: these add to another existing observation performed in 2006 
(OBSID: 07745, PI: Bianchi) to reach an exposure time of ~ 110 ks, allowing us to study specific 
features of the soft X-ray emission (0.3 - 2 keV) from this source. 

There is a remarkable correspondence between the soft X-ray and the bi-conical NLR and 
ENLR as mapped by [O m] emission, and the radio jets, as originally reported by Bianchi et 
al. (2010) and Gonzalez-Martin et al. (2010), suggesting a physical interplay among the three 
components. This is even more evident in the comparison of the merged image of all available 
Chandra observations with the HST and VLA images presented in Figure 1 . 

We have performed spectral analyses on different region of the source, and modeled the ex- 
tracted spectra with different fitting models. We find that the nuclear emission is best fitted with 
a highly photoionized plasma (log U = 1.1+0.1) and shows a significant (4<x) decrease of the 
observed soft X-ray flux that can be interpreted as the result of intrinsic variability of the nuclear 
continuum. The extended cone emission is effectively modeled with a mildly photoionized plasma 
(log U = 0.9 + 0.3 for the NW cone and log U = 0.7 + 0.2 for the SE cone) diffused over the bicone. 
A collisionally ionized gas at a temperature ~ 0.8 keV appears to be located ahead of the observed 
radio knots, at the interface between the optical arcs. Another denser, hotter thermal component 
at ~ 1.1 keV is located between the nucleus and the radio knots. The thermal components can 
be interpreted in terms of shock interaction of the radio ejecta with the galaxy plane, with the jet 
depositing > 30% of its power into the ISM. This represents a > 0.05% fraction of the available 
accretion power that is thermally coupled to the ISM, substantially lower than what required by 
the majority of quasar feedback models, but compatible with a two-stage feedback. 

We detect diffuse, soft X-ray emission extending up to ~ 13 kpc along the bicone axis. In the 
region perpendicular to this axis, where the obscuring torus is supposed to prevent nuclear X-ray 
reaching the ISM, we find significant emission extending up to ~ 9" from the nucleus. We model 
this emission with a plasma featuring a lower ionization parameter with respect to the emission 
extracted in the bicone region, an so we interpret the former as originating from the nuclear contin- 
uum filtered by a warm absorber with an uncovered factor < 3%. Deeper observations are however 
required in order to better constrain the model parameter in this region of faint emission. 
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Fig. 1. — (Upper- left panel) PSF-deconvolved ACIS image (0.3 - 2 keV) using the EMC2 method, 
with subpixel binning (1/8 of the native pixel size) and a 2X2 FWHM gaussian filter smoothing. 
The same image is presented in the lower-left and lower-right panels with overlaid the HST [O m] 
and radio VLA 6 cm contours, respectively. (Upper-right panel) HST [O m] image. 
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Fig. 2. — Radial profiles of the ionizing cones (black dots) in comparison with the simulated PSF 
(red dots) in the 0.3 - 2 keV range. 
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A. Spectral fits results 

In this Appendix we present results of spectral fits for the models presented in Sects. 4.1 
and 5.2, both for the cone (Tables 4, 5 and Figures 12, 13) and the nuclear (Tables 6, 7, 8 and 
Figures 14, 15) regions. When possible, fits were performed each observation as well as on merged 
observations. 
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Table 1 : Physical parameters of the collisionally ionized gas found 
in regions of increased ionization represented in Figure 8. 



Region 


1 


2 


3+4 


n e (cm -3 ) 


0.48 + 0.10 


0.52 + 0.12 


0.14 + 0.02 


Pth (10~ 9 dyne cm -2 ) 


1.78 + 0.41 


2.34 + 0.61 


0.38 + 0.08 


E th (10 53 erg) 


0.79 + 0.18 


1.04 + 0.26 


2.73 + 0.57 


tcooi (10 7 yr) 


0.76 + 0.32 


1.24 + 0.64 


1.63 + 0.53 


c s (km s _1 ) 


668 + 103 


737 + 128 


565 + 73 


v sh (km s _1 ) 


905 + 46 


999 + 59 


765 + 53 


t cross (10 5 yr) 


3.66 + 0.57 


3.32 + 0.58 


4.33 + 0.56 


E t h It cro SS (10 41 erg s" 1 ) 


0.69 + 0.19 


1.00 + 0.30 


2.00 + 0.49 



Table 2: Observations properties 



Obs. ID 


07745 


12294 


13124 


13125 


merged 


Obs. Date 


2006-11-18 


2010-09-16 


2010-09-17 


2010-09-19 


- 


Exposure (ks) 


35.07 


9.92 


52.37 


16.76 


114.12 


Net Counts 0.3 - 2 keV (error) a 


3173(57) 


739(27) 


3649(61) 


1216(35) 


8777(95) 


Net Counts 2-10 keV (error) a 


345(23) 


79(12) 


443(27) 


160(15) 


1027(41) 



Notes: 

a Counts evaluated in a circular region of 15" radius. 
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Fig. 3. — Adaptive smoothing applied to subpixel binned ACIS-S image of the soft X-ray emission 
(0.3 - 2 keV); showing extended emission up to ~ 6 kpc. Contours correspond to seven logarithmic 
intervals in the range of 0.003 - 5% with respect to the brightest pixel. 



Table 3: Best fit photoionization models for other regions 


(merged observations). 




Region 


1 


2 


3+4 


[O m] outer arcs 


cross-cone region 


Net Counts 0.3-2 keV (error) 


1002(32) 


1248(36) 


691(26) 


308(18) 


547(24) 


Model Parameter 












logtfi 


0.87*$ 


0.85 ±0.12 


1.33 ±0.31 


0.183g 


n o i +0.08 
U - 81 -0.05 


log N H i 


21.41!™ 


21.55 ±0.58 


20* 


20.07!gj5 


20* 


F\ (0.3-2) a 


0.22 ± 0.06 


19 +004 
ua ^-0.03 


0.04 ± 0.03 


0.17 ±0.02 


0.16 ±0.03 


\ogU 2 


i 1 7+O.6O 
11 '-0.33 


-0.75 ± 0.28 


-0.85!^ 


-2* 


< -2.42 


\ogN H 2 


20* 


21* 


20* 


20* 


20* 


^2 (0.3-2) " 


0.13 ±0.07 


16 +004 
u>1 °-0.03 


26 +003 

u - ZD -0.06 


0.06 ± 0.03 


04 +003 
u - u -0.02 


kT b 


1.06_ 011 


1.29 ±0.15 


0.76 ±0.11 


- 


- 


Fc (0.3-2) a 


0.05 ± 0.02 


04 +002 
u - u ^-o.oi 


0.08!^ 


- 


- 


* 2 (dof) 


0.98(20) 


0.59(20) 


0.72(20) 


0.98(9) 


1.01(20) 


-Fo.3-2 a 


0.39 ± 0.04 


0.39 ± 0.02 


0.42 ± 0.03 


0.23 ± 0.02 


29 003 

"•^^-0.02 



Notes: 

a Unabsorbed flux in units of 10 

b Plasma temperature in keV. 



erg cm 2 s l . 
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Fig. 4. — Extended radial profiles of the NW (top panel) and SE (bottom panel) cones (black dots) 
in comparison with the simulated PSF (red dots) in the 0.3 - 2 keV range. 
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Fig. 5. — Adaptive smoothing applied to subpixel binned ACIS-S image of the soft X-ray emission 
(0.3 - 2 keV) with overlaid the extraction regions considered in the text. 
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Fig. 6. — Radial profiles of the cross cone region (black dots) in comparison with the simulated 
PSF (red dots) in the 0.3-2 keV range. 
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Fig. 7. — Emission-line maps of the central 3.5 kpc of Mrk 573: from left to right are shown the O 
vii triplet (0.53 - 0.63 keV), O vm Lya (0.63 -0.68 keV) and Ne ix triplet (0.90 - 0.95 keV) maps, 
with overlaid the 6 cm radio contours (upper panel) and [O m] contours (lower panel). The images 
are rebinned to 1/8 of the native pixel size and smoothed with a 3X3 FWHM gaussian filter. The 
nucleus position is marked with the black cross. 
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Fig. 8. — Ratio Ne ix / O vn maps (left panel) and Ne ix / O vm (right panel) with overlaid [O m] 
(subpixel binning 1/2 of the native pixel) smoothed with a 3X3 FWHM gaussian filter. The high 
ionization regions used for spectral extraction are shown with white full lines. The with dashed 
circle represents the 1" radius region used for the nuclear spectrum extraction, while the blue 
dashed circle represents the inner 0.5" radius region used for restricted nuclear fitting discussed in 
Sect 4.1. 
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Fig. 9. — Variation of the [O m] to soft X-ray flux ratio, in the NW (upper panel) and SE (lower 
panel) cone as a function of the distance from the nucleus. The vertical dashed lines represent 
the location of the inner and outer optical arcs, and the squares represent the regions of increased 
ionization shown in Figure 8. 
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Fig. 10. — Merged Chandra/ 'ACIS-S 0.3 - 2 keV spectra of the regions discussed in Sect. 4.2 
and 5, with best fit photoionization model; for each spectrum we also show with dotted lines 
the additive components of the best fitting models. From top-left to bottom-right we show fits 
to spectra extracted in region 1, 2 and 3+4 shown in Figure 8, regions of outer optical arcs, and 
cross-cone region shown in Figure 5. 
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Fig. 11. — Variations of the nuclear spectral components between 2006 obs. 07745 and 2010 
CHEERS obs. Black points represent emission line fluxes (triangles indicate upper limits), the 
blue point represents the power-law normalization, the red point represents the soft X-ray flux and 
the green point represents the hard X-ray flux. Emission line fluxes and the power-law normal- 
ization are expressed in units of 10~ 5 photons cm" 2 s" 1 and X-ray fluxes are expressed in units of 
10 -13 erg cm -2 s _1 . 
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Table 4: Measured line fluxes of the cone regions 



Obs. ID 




07745 


13124 


CHEERS OBS 


ALL OBS 


Cone 


NW 


SE 


NW 


SE 


NW 


SE 


NW 


SE 


Net Counts 0.3 


2 keV (error) 


358(19) 


277(17) 


424(21) 


371(19) 


656(26) 


550(24) 


1013(32) 


826(29) 


Line 


Energy 11 


Flux" 


CVHey 


0.371 




- 














N VI triplet 
C VI Ly/3 


0.426 
0.436 


q °6+ - 62 
0.40+8$ 

165 +8:1 


2.03 + 0.46 


0.75 ± 0.36 


1.73+0.35 


0.71+0.36 


1.60 + 0.32 


0.74 + 0.31 


1.75 + 0.26 


N VII Lya 


0.500 


- 


<0.48 




<0.46 




0.30 + 0.20 




VII triplet 


0.569 


1.65 + 0.27 


1.16 + 0.25 


1.42 + 0.21 


1.14 + 0.24 


1.45+0.19 


1.30 + 0.20 


1.52 + 0.15 


VIII Lya 


0.654 


<0.32 


0.26 + 0.13 


0.38 + 0.11 


0.40 + 0.11 


0.38 + 0.10 


0.44 + 0.10 


0.32 + 0.08 


0.38 ± 0.08 


O VII Hey 


0.698 




- 














Fe XVII 3s2 


0.727 




- 














VII RRC 
Fe XVII 3d2p 


0.739 
0.826 


31+ 014 
" -0.12 

<0.21 


0.40 + 0.11 
0.23 ± 0.07 


0.43 ± 0.09 
0.18 + 0.11 


0.27+0.08 
0.11+0.05 


48+ 009 

-848 

20+°* 


0.26 ± 0.08 
0.09 ± 0.05 


0.42 ± 0.07 
0.17+0.09 


0.31 ±0.06 
0.13 + 0.04 


VIII RRC 


0.871 


43+0.11 
"'+ -0.25 


- 


<0.19 




< 0.15 




<0.23 




Ne IX triplet 
Fe XX 3<Cp 


0.915 
0.965 


0.13 + 0.06 


22+ 007 
"■ zz -0.10 


0.18+0.05 


24+ 005 
"■ z +-0.09 


0.17+0.04 


0.23+0.08 


0.16 + 0.04 


Ne X Ly<? 


1.022 


0.08+j™ 

06+ ii ™ 
"" -0.05 


0.09 ± 0.05 


0.13 + 0.03 


0.06 ± 0.03 


0.14 + 0.03 


0.07 ± 0.03 


0.14 + 0.03 


0.07 ± 0.03 


Ne IX Hey 
Ne IX He* 


1.127 
1.152 


- 


0.04 ± 0.02 


<0.03 


0.04 ± 0.02 


<0.03 


0.05 ± 0.02 


<0.03 


NeXLyS 


1.211 




- 














Mg XI triplet 


1.352 


<0.08 


- 


0.06 ± 0.03 




0.05 ± 0.03 




0.05 ± 0.02 




Fe XXII 4p2p 


1.425 




- 














Si XIII triplet 1.839 
Power-Law norm. 


66+ 016 

"™-n 55 


0.31+0.08 


0.31+0.09 


44+0.07 
"■++-0.08 


0.30 + 0.11 


0.43 ± 0.07 


0.32 + 0.11 


0.37 ± 0.05 


^(dof) 
^0.3-2° 


0.72(3) 

n ro+0.04 

,I58 -0.06 


0.82(4) 
0.51+0.04 


1.19(6) 
0.45 ± 0.03 


0.77(8) 
0.49 ± 0.03 


0.88(16) 
0.45 ± 0.03 


0.62(16) 

0.48+0.03 


1.29(31) 
0.47+0.02 


0.75(28) 
0.49 ± 0.02 



Notes: 
a Line rest-frame energy in keV. 

Line fluxes in units of lO -3 photons cm s . 
c Unabsorbed flux in the 0.3 - 2 keV band in units of 10~ 13 erg cm -2 s _1 . 



Table 5: Best fit photoionization models for the cone regions 



Obs. ID 


07745 


13124 


CHEERS OBS 


ALL OBS 


Model Parameter 


NW 


SE 


NW 


SE 


NW 


SE 


NW 


SE 


log Ui 


1.00 + 0.08 


75 +UJ0 
U -' J -0.19 


0.78 + 0.15 


6 2+0.29 
"• D -0.23 


75+ 0J4 
"■"-0.19 


55+0.23 
""'-'-0.18 


0.93 ± 0.28 


69+ IIJI 
"■"'-0.19 


loglVff! 


20* 


20* 


20* 


20* 


20* 


20* 


20* 


20* 


F t (0.3-2) a 


36+ 004 
" -0.05 


26+ 008 


025 +0.13 

u -0.0? 

-0 55+0-34 

"••'-'-0.39 


2 7+0.09 
"■ z '-0.11 


024+0.11 
"■ z +-0.07 


27+0.09 
"■ -0.08 


21+0- 10 
"■ zl -0.07 


0.22 ± 0.08 


log(/ 2 


-0.5* 


-1* 


-0 95+ - 32 
"■' J -0.26 


-0 50+ 055 
" J "-0.33 


n 05+0.40 
-0.85_ 027 


_ 051 +0.27 
"""-0.31 


-0.80+^1 


log« H2 


20* 


20* 


20* 


20* 


20* 


20* 


20* 


20* 


f2 (0.3-2)" 


32+ 006 


0.41 ±0.11 


26+012 
" zo -0.08 


34+0.12 
"■ -0.11 


025+011 
"■ -0.07 


30+0.09 
"■ -0.10 


31+ 009 
"■ J1 -0.07 


0.37+0.08 


kT h 






0.82 + 0.13 


1 04+0- 23 
'■"+-0.17 


0.89 + 0.11 


1 09+0- 41 


0.80 + 0.09 


0.97+0.17 


^C(0.3-2) a 






07 +0.05 
""'-0.04 


0.04 ± 0.03 


07+ 003 
"■"'-0.05 


0.04 ± 0.02 


O.08_ o ; 03 


0.04 ± 0.03 


^(dof) 


1.16(12) 


0.96(9) 


1.04(12) 


0.47(11) 


1.10(22) 


0.42(19) 


1.20(37) 


0.57(31) 


Fo.1-2" 


n 6S+ 004 
0.68_ nm 


0.68 ± 0.05 


57 +0.05 
"■^'-O.M 


64+006 
"■"+-0 05 


0.57 ± 0.03 


61+004 


0.60 ± 0.04 


0.62 ± 0.04 



Notes: 

a Unabsorbed flux in the 0.3 - 2 keV band in units of 10~ 
Plasma temperature in keV. 
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Fig. 12.— Chandra/ ACIS-S 0.3 - 2 keV spectra of the NW (left column) and SE (right column) 
cones emission (extracted from the region shown in Figure 5) in different observations, with best fit 
emission lines models. From top to bottom we show fits to OBS. 07745, 13124, merged data from 
CHEERS observations and merged data from all observations. For spectra of single observations 
we also show with dotted lines the additive components of the best fitting models. 
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Fig. 13.— Chandra/ ACIS-S 0.3 - 2 keV spectra of the NW (left column) and SE (right column) 
cones emission (extracted from the region shown in Figure 5) in different observations, with best fit 
photoionization models. From top to bottom we show fits to OBS. 07745, 13124, merged data from 
CHEERS observations and merged data from all observations. For spectra of single observations 
we also show with dotted lines the additive components of the best fitting models. 
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Table 6: Measured line fluxes of the nuclear region 



Obs. ID. 




07745 


12294 


13124 


13125 


CHEERS OBS 


ALL OBS 


Net Counts 0.3 


- 2 keV (error) 


2109(46) 


460(21) 


2308(48) 


765(28) 


3533(59) 


5642(75) 


Line 


Energy a 


Flux b 


CVHey 


0.371 


1 C -7T + 10.23 
18./3_ 740 


- 


12.22! 2 ; 2 5 


- 


11.88* 212 


19710.55 


N VI triplet 


0.426 


- 


- 


- 


- 


- 


- 


C VI Ly/3 


0.436 


<6.84 


7 S7+ 2 - 51 
'••"-2.31 


o en+0-86 
Z - 8U -0.95 


6.18!&» 


9 07+O.8I 
z -°'-0.83 


-j «+0.18 
3 - 65 -0.23 


N VII Lya 


0.500 


5 33+ 111 

J '-'- , -2.08 


<2.26 


2.29 ± 0.53 


4.10!- 9 


2.67 ± 0.45 


-1 9O+0.13 
3.Z8_ 012 


O VII triplet 


0.569 


9 69 +107 


7 02+ L69 

'■ uz -1.38 


6.90 ±0.61 


7 40 +L47 

' - w -4.03 


7.00 ± 0.50 


7.76 ±0.13 


O VIE Lya 


0.654 


3.14 ±0.43 


1 74+ 113 


i c/c+0.32 
L56 -0.28 


2.47 ± 0.60 


1 70+0.25 
1 ■ ' ^-0.22 


9 . 4 +0.05 
^ -0.07 


O VII Hey 


0.698 


- 


f\ *; 1+0-47 
°- 61 -0.54 


- 


- 


- 


- 


Fe XVII 3s2 


0.727 


- 


- 


<0.43 


- 


- 


- 


O VII RRC 


0.739 


2.09 ± 0.34 


- 


1 34+0.43 

i - r, ^-0.21 


0.55 ± 0.41 


1 21 +017 

i ' z,1 -0.22 


1 4C+0.03 

1/+J -0.06 


Fe XVII 3d2p 


0.826 


1 Q/c+0-36 
L86 -0.50 


1 36 +039 


1 10 +017 

iau -0.20 


1.34 ±-0.36 


1 1A+0.19 
116 0.14 


1 oT+0.04 
i -°-'-0.03 


O VIII RRC 


0.871 


< 1.02 


o.8o:^ 3 8 


- 


- 


<0.19 


11+001 


Ne IX triplet 


0.915 


9 94+0.44 
^■^ -0.62 


1.083£ 


i /ci +0.18 
161 0.15 


i 4Q+0.36 
1 ,H ^-0.44 


1 CO +0.22 

1.38_ 013 


1 7 j; +0.04 
i - /J -0.03 


Fe XX 3d2p 


0.965 


<0.60 


1 oq+0.56 
L/8_ 03 4 


23 +014 


<0.53 


39+O.II 
U - JZ -0.19 


32 +002 

""^ -0.03 


Ne X Lya 


1.022 


1.10 ±0.20 


<0.32 


35+0- 10 


0.82! ; 29 


39+ 012 


55 +002 

u - JJ -0.03 


Ne IX Hey 


1.127 


0.44 ±0.12 


41 +0 ' 22 


22 +010 

u - zz -0.08 


35 +018 

U - JJ -0.17 


33+ 012 

u -"-0.07 


0.40 ± 0.02 


Ne IX He6 


1.152 


- 


- 


031 +0.15 


- 


1 5 +a2 ° 

U ' iJ -0.07 


0.07 ± 0.02 


Ne X LyyS 


1.211 


0.51 ±0.11 


041 +o.i8 
u/+1 -o.io 


27 +007 

"•^'-0.08 


0.28 ±0.15 


30 +009 


038!$ 


Mg XI triplet 


1.352 


0.53 ±0.10 


0.59 ±0.18 


0.49 ± 0.07 


0.23 ±0.13 


0.44 ± 0.06 


0A71Z 


Fe XXII 4p2p 


1.425 


0.16 ±0.09 


0.18*$ 


(MB** 


- 


n n^+0.05 
0.06_ 004 


09 +002 
u - uy -o.oi 


Si XIII triplet 


1.839 


0.32 ± 0.08 


- 


0.21 ± 0.05 


0.15 ±0.12 


0.19 ±0.04 


0.23 ± 0.01 


Power-Law norm. 


3.68 ±0.51 


4.48_^ 


4 20 +032 

H ' z ' u -0.30 


4.71 ±0.74 


4 28 +0 ' 25 


4 - 06 -0.06 


* 2 (dof) 




0.93(49) 


1.28(4) 


1.28(55) 


0.54(13) 


1.21(102) 


1.63(169) 


-F(0.3-2) C 




5 40 +0 - 29 

J - w -0.23 


3.69 ± 0.22 


4.01 ±0.15 


3 7Q+0.82 


40 c+0.13 


4.37 ±0.11 



Notes: 

a Line rest-frame energy in keV 

b Line fluxes in units of 10~ 5 photons cirT 2 s~ 



c Unabsorbed flux in the 0.3 



2 keV band in units of 10 13 erg cm 2 s '. 
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Fig. 14. — Chandra/ ACIS-S 0.3 - 2 keV spectra of the nuclear emission (extracted from the region 
shown in Figure 5) in different observations, with best fit emission lines models. From top left 
to bottom right we show fits to OBS. 07745, 12294, 13124, 13125, merged data from CHEERS 
observations and merged data from all observations. For spectra of single observations we also 
show with dotted lines the additive components of the best fitting models. 
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Table 7: Best fit photoionization models for the nuclear region 



Obs. ID. 

Net Counts 0.3 - 10 keV (error) 


07745 
2352(48) 


12294 
508(23) 


13124 
2663(52) 


13125 
840(29) 


CHEERS OBS 

4063(64) 


ALL OBS 
6415(80) 


Model Parameter 














logt/i 

logN/n 

F\ (0.3-2)" 1 
•Fl(2-10) d 

\ogU 2 

log/V H2 

^2(0.3-2) d 
^2(2-10) d 

Refl c 

FsiKa" 

EsKa 

FsKa B 

EFeKa 

FFeKa" 


20.62!°i° 

9 77+0.19 
^•"-0.18 
041+ 0.05 
u - H1 -0.05 

-0 75 +013 
20.22!°;g 
2.86! ; 28 

24 +007 

u ' ZH -0.05 

1 20 +016 
1.839* 

0.14 ±0.07 

2.430* 

0.18 ±0.08 

6.35 ±0.01 

0-89^ 


i i i +0.09 
1,1 -0.07 

21 20 +102 

9 7^+0.40 
z - /J -0.45 
A A/1 +0.33 

0.64_ 021 
-0.89 ± 0.63 
21.86 ±1.60 

i 17 +0.56 
i - i/ -0.53 

34 +0 - 74 

"•-' -0.22 

70 +0 - 37 

U -' U -0.51 

1.839* 

2.430* 
<0.41 


l.OS ^ 
20.99!"; 34 

1 Q4+0.12 

0.30 ± 0.05 

-0 73 +a28 

20.82+ ; 38 
i 47+0.22 

A 1^+007 

O- 1 6 -o.04 
0.69 ±0.19 

1.839* 
0.10 ±0.05 

2.430* 
0.13 ±0.06 

£ ^o+0.09 
°- J:, -0.06 

1 20 +0 - 26 

1 - zu -0.23 


1 09 +009 
20 91 +0 ' 44 

zu - yi -0.32 

2.19 ±0.28 
37+ 010 

u --"-0.09 

-0.88 ± 0.23 

21.66!?- 73 

2.10t| 

G.48^ 3 

0.73 ± 0.65 

1.839* 

<0.22 

2.430* 

<0.22 

A m+100 

6.03_ 003 
9V - 52 

U ' y:, -0.43 


1 05 +004 

i - UJ -0.03 

20.99!^| 

2 03 +oi2 

fl 39+0.04 
u - JZ -0.05 

-0 72 +0 ' 26 

21.04 ±0.60 

1 51 +0 ' 21 

1,J1 -0.20 

°- 18 -0.03 

0.69!™ 

1.839* 

0.08 ± 0.05 

2.430* 
0.11 ±0.05 

^ co+0.10 
°" :, - , -0.07 

1 13 +0 ' 24 

laJ -0.20 


1.06 ±0.03 
20.88!™ 

9 91 +0.14 
^•^ -0.13 

34 +005 

u - J ^-0.04 
n 77+0.13 
v.i i_ ou 

20.70!" 34 

1 86+ - 20 

20 +006 

u - zu -0.04 

0.82 ±0.14 
1.839* 

0.10 ±0.04 
2.430* 

0.13 ±0.04 

A /I +0.06 
6 - 48 -0.05 

1 14+0.19 

i - iH -0.17 


* 2 (dof) 

-F(0.3-2) d 

-F(2-io) a 


0.94(68) 

-r 77+0.22 
J -''-0.21 

3.84 ± 0.32 


0.79(14) 
4 00 +a29 
2.36_ 057 


1.42(76) 

4.06 ±0.14 

3 64 +031 

- , -°^-0.29 


1.10(23) 
4 40 +0 - 35 

o 09+0.46 
- , - ,z -0.60 


1.21(132) 
4.09 ±0.13 

o -to +0.28 
- , "-'- , -0.25 


1.55(211) 
4^4+0.11 

H - J -0.06 

3.68 ± 0.21 



Notes: 

a Unabsorbed flux in units of 10~ 13 erg cirT 2 s~' . 

b Plasma temperature in keV. 

c Normalization of the reflection component in units of 10~ 3 photons keV -1 cm -2 s -1 



Line rest-frame energy in keV. 
e Line fluxes in units of 10 photons crrT 2 s~ 
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Fig. 15. — Chandra/AClS-S 0.3 - 10 keV spectra of the nuclear emission (extracted from the 
region shown in Figure 5) in different observations, with best fit photoionization models. From top 
left to bottom right we show fits to OBS. 07745, 12294, 13124, 13125, merged data from CHEERS 
observations and merged data from all observations. 
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Table 8: Best fit photoionization models for the 


inner nuclear region 


Obs. ID. 


07745 


13124 


CHEERS OBS 


Net Counts 0.3 - 10 keV (error) 


1307(36) 


1441(38) 


2030(45) 


Model Parameter 








logC/i 


1.10 + 0.06 


1.02 + 0.05 


1 05 +004 

1 - UJ -0.05 


logNtfi 


20 51 +032 


20 95+ 021 


on 99+ 016 
z,u.^_ 029 


^1 (0.3-2) 3 


1.55 + 0.17 


1.10 + 0.09 


1 12 +(m 


Fl(2-W) a 


r> 9 9 +0.06 


n 1 £+0.04 

°- 16 -0.03 


1 7 +0 - 02 

ua/ -0.04 


logU 2 


-°- 68 -o:24 


-0 79 +0 - 39 


n '7^+°- 31 
-0.76_ 017 


logN H2 


20 10 +0 - 48 


20 77 +0 - 52 

z ' u - / '-0.43 


21.06 + 0.86 


^2(0.3-2) a 


1.61+0.21 


77 +018 

V. I /_ 01 4 


0.80!- 3 


^2(2-10) a 


12 +005 


09 +004 


o.ioSS 


Refl c 


0.56_ 019 


40 +013 


40 +012 

u - HU -0.07 


F d 


1.839* 


1.839* 


1.839* 


FsiKa S 


0.11+0.05 


0.10 + 0.04 


0.07 + 0.03 


T7 d 


2.430* 


2.430* 


2.430* 


FsKa £ 


0.08 + 0.07 


0.06 + 0.04 


0.06 + 0.04 


F d 


-r O-7+0.03 
6 - 37 -0.02 


r ro+0.08 

O.D8_ 007 


6.58 + 0.07 


FFeKa £ 


0.60t\ 6 8 


0.78SS 


0.80SS 


^ 2 (dof) 


0.98(47) 


1.26(57) 


1.19(95) 


^(0.3-2) a 


o 0^+0.17 

3.36_ 012 


2.43 + 0.09 


2 14 +0 - 07 


^(2-10) 3 


2.21 + 0.25 


i ao+0.15 

Z - J8 -0.19 


9 -3 1+0.18 
z " : ' 1 -0.10 



Notes: 

a Unabsorbed flux in units of 10" erg cm 2 s -1 . 

b Plasma temperature in keV. 

c Normalization of the reflection 

10~ 3 photons keV" 1 cm 2 s _1 

d Line rest-frame energy in keV. 

e Line fluxes in units of 10~ 5 photons cm 2 s 



component in units of 



